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ABSTRACT 

We have used the IRAM Plateau de Bure millimeter interferometer and the UKIRT 1- 
\Q . 5 /im Imager Spectrometer (UIST) to test the connection between the major phases of 

' spheroid growth and nuclear accretion by mapping CO emission in nine submillimetre- 

detected QSOs at z = 1.7-2.6 with black hole (BH) masses derived from near-infrared 
■ spectroscopy. When combined with one QSO obtained from the literature, we present 

OO ' sensitive CO(3-2) or CO(2-l) observations of 10 submillimetre-detected QSOs se- 

. lected at the epoch of peak activity in both QSOs and submillimetre (submm) galaxies 

(SMGs). CO is detected in 5/6 very optically luminous (Mb ~ —28) submm-detected 
QSOs with BH masses Mbh — 10 9 -10 10 M Q , confirming the presence of large gas reser- 
voirs of M gas ~ 3.4 x 10 10 Mq. Our BH masses and dynamical mass constraints on 



$-i ' the host spheroids suggest, at face value, that these optically luminous QSOs at z = 2 

lie about an order of magnitude above the local BH-sphcroid relation, AfBH/M sp h, al- 
though this result is dependent on the size and inclination of the CO-emitting region. 
However, we find that their BH masses are ~ 30 times too large and their surface 
density is ~ 300 times too small to be related to typical SMGs in an evolutionary 
sequence. Conversely, we measure weaker CO emission in four fainter (Mb ~ —25) 
submm-detected QSOs with properties, BH masses (Mbh — 5 x 10 8 M q ), and surface 
densities similar to SMGs. These QSOs appear to lie near the local MBH/M sp h rela- 
tion, making them plausible 'transition objects' in the proposed evolutionary sequence 
linking QSOs to the formation of massive young galaxies and BHs at high-redshift. 
We show that SMGs have a higher incidence of bimodal CO line profiles than seen 
in our QSO sample, which we interpret as an effect of their relative inclinations, with 
the QSOs seen more face-on. Finally, we find that the gas masses of the four fainter 
submm-detected QSOs imply that their star formation episodes could be sustained 
for ~ lOMyr, and are consistent with representing a phase in the formation of mas- 
sive galaxies which overlaps a preceding SMG starburst phase, before subsequently 
evolving into a population of present-day massive ellipticals. 

Key words: galaxies: high-redshift - galaxies: evolution - galaxies: formation - 
galaxies: kinematics and dynamics - quasars: emission lines - submillimetre 
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1 INTRODUCTION 

It has been established that every massive, local spheroid 
harbours a supermassive black hole (SMBH) in its cen- 
tre whose mass is proportional to that o f its host (e.g. 
iMagorrian et al.lll99sl ; IGebhardt et alj|2000l ). This suggests 
that the black holes (BHs) and their surrounding galaxies 
were formed synchronously. This hypothesis has found sup- 
port from hydrodynamical simulations of galaxy formation, 
which use feedback from winds and outflows from active 
galactic nuclei (AG N) to link the growth of the SMBH to 
that of its host (e.g. |Pi Matteo. Springel fc Hernquistll2005l ; 
iHopkins et"aD 120051 ; iBower et al.ll2006l). Thus these model s 
support a picture, first presented by ISanders et al] <|l988h . 
where a starburst-dominated ultra-luminous infrared galaxy 
(ULIRG), arising from a merger, evolves first into an ob- 
scured QSO and then into an unobscured QSO, before fi- 
nally becoming a passive spheroid. 

The high-redshift population of ULIRGs in this 
proposed evolutionary cyc le are the submil li metre 
(submm) galaxies (SMGs; ISmail. Ivison fc Blair] 1 19971 ; 
IChapman etail 120051 ; ICoppin et alj|2006l ). These systems 
have ULIRG-like bo l ometric luminositie s, Lm > 10 12 L© 
l|Kovacs et al] 120061 ; ICoppin et al] 120081 ). and they have 
many of the p r opert ie s expected for gas-rich m ergers 
i|Swinbank et al] 12004 120061 ; iTacconi et ail |2006h . This 
population evolves rapidly out to a p eak at z ~ 2.3, crudely 
matching the evolution of QSOs l|Chapman et al] 120051 ) 
and providing additional circumstantial evidence for a 
link between SMBH growth and spheroids. Two further 
results have shed light on the evolutionary link between 
SMGs and QSOs. Firstly, a modest fraction of optically 
luminous QSO s at z ~ 2 are dete cted in the submm/mm 
(~ (25 ± 10)%; IOmont et alj|2003h showing that the QSO- 
and SMG-phases do not overlap significantly, given the 
lifetime estimates of the two populations (QSOs make up 
~ 4% of flux-limited samples of SMGs; IChapman et al] 
2005). But when a QSO is detected in the submm/mm then 
it could to be in the transition phase from an SMG to an 
unobscured QSO, making it s properties a powerful probe of 
the e volutionary cycle (e.g. IStevens et al] 12005; Pag e et al] 
12004! ). Secondly, the evolutionary state of the SMBHs 
within SMGs can also be judge d using the 2-Ms Chan dra 
Deep Field North observations (|Alexander et al] l2003h to 
derive accurate AGN luminosities and hence lower limits on 
the BH masses (Mbh) in those SMGs wi th precise redshift s 
in this region (|Alexander et al.ll2005aH bl; iBorvs et al]|2005h . 
These studies suggest that the AGN in typical SMGs are 
growing almost continuously - but that the SMBHs in these 
galaxies appear to be several times less m assive than seen 
in co mparably massive galaxies at z ~ (|Alexander et al] 
l2008h . 

Together these results argue for a fast transition 
from a star- formation-d ominated SMG-p hase to the AGN- 
dominated QSO-phase (|Page et al.ll2004r ). The latter phase 
will result in the rapid BH growth necessary to account 
for the pres ent-day relation between spheroid and SMBH 
masses (e.g. IMagorrian et al.l fl99Sl ; IGebhardt et al] l200d) . 
Can we confirm this and more generally test the proposed 
evolutionary link between SMGs and QSOs at the peak of 
their activity at z ~ 2 ? 

This evolutionary cycle has been tested in the local Uni- 



verse bv_romrjaxmgthe properties of QSOs and ULIRGs 
(e.g. ITacconi et al]|2002l ). However, both of these popula- 
tions are 1000 x less abundant in the local Universe than 
they were at the era of their peak activity at z ~ 2 and 
so we have to be cautious about extrapolating from local 
examples to the high-re dshift progenitors of the bulk of to- 
day's massive spheroids (|Genzel et al.| [2003 ; Swi nbank et al.l 
2006). Thus, to properly test the validity of this cycle for 
typical spheroids we have to compare QSOs and ULIRGs at 
the era where their populations peaked: z ~ 2. 

The critical pieces of information needed to test the 
link between SMGs and QSOs are the relative dynamical, 
gas and SMBH masses of these two populations. In prin- 
ciple the dynamical masses can be derived from optical or 
near-infrared observations of emission line gas in the SMGs 
or QSOs (see ISwinbank et alJliooi . 120051 . I2OO6I ). However, 
there is the problem of removing the QSO emission in these 
observations, as well as the effects of dust obscuration and 
outflows. In contrast, molecular CO emission line profiles 
are relatively immune to the effects of obscuration and out- 
flows, while at the same time yielding additional constraints 
on the relationship between QSOs and SMGs from their gas 
masses. 

There is currently a lack of sensitive CO observations of 
QS Os at z ~ 1-3, with data published on only eight sources 
(e.g. lFraver et al.lll99Sl; iGuilloteau et alj|l999 : iBeelen et al] 
|2004 lHainline et alj|2004 . Instead the focus has been on 
CO studies of QSOs at z ;> 4 (e.g. lOmont et al] Il996l ; 
IWalter et ail |2004 iRiechers et~aTl 120061 ). although these 
QSOs have little overlap with the redshift range where SMGs 
are typically detected. The paucity of CO constraints for 
z > 1 QSOs also reflects the difficulty in determining their 
systemic redshifts with sufficient precision to guarantee that 
the CO emission falls within the bandwidth of typical mil- 
limetre (mm) correlators. However, sensitive near-infrared 
spectroscopy of the Civ, Mgn, and [Oni]5007 emission lines 
in QSOs can provide redshifts with required precision as well 
as Ha or H/3 fluxes and line widths to yiel d Mbh estimates 
(|Takata et al]|2006l ; [Alexander et alj|2008l) . 

We have carried out a quantitative test of the proposed 
link between SMGs and QSOs at 2 ~ 2 where both popula- 
tions are most common. We have obtained precise systemic 
redshifts from near-infrared spectroscopy of potential transi- 
tion QSOs (i.e. submm/mm-detected QSOs) and then used 
the IRAM Plateau de Bure Interferometer (PdBI) to search 
for CO emission. We relate their dynamical, ga s and SMBH 
mass es to SMGs from the PdBI CO survey (|Greve et al] 
2005). We test: a) whether the cold gas masses in these QSOs 
are similar to those in SMGs; b) whether the line widths and 
dynamical masses of these two populations are comparable; 
and c) how the ratio of SMBH to dynamical masses for these 
submm-detected QSOs relate to the estimates for SMGs and 
those for optically lumin ous QSOs (which lie on t he present- 
day Mbh-M ct relation; iMcLure fc Dunlod I200I ) . Together 
these observations can constrain the proposed evolution- 
ary sequence which links QSOs to the formation of massive 
young galaxies and SMBHs at high redshift. 

We describe the sample selection, observations and data 
reduction in 3U The results of the near-infrared and mm 
CO spectra are given in 33 The CO properties of the 
submm-detected QSOs are compared and contrasted with 
SMGs in 331 and we discuss the evolutionary status of the 
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Figure 1. R-band apparent magnitude versus redshift for our 
sample of submm-detected QSOs. We show that these lie within 
the interquartile range (shaded region) of the redshift distribution 
of SM Gs (median <z>~ 2.3; dashed line) from [Chapman et al.l 
(2005). This shows that we are probing a wide range in optical 
luminosities (and potentially a wide range in A/bh) at the epoch 
where the QSO and SMG populations peak. 

submm-detected QSOs in ij5] Our conclusions are given in 
Sj6] W e adopt cosm o logica l parameters from the WMAP 
fits in ISpereel et all (|2003l ): ft A = 0.73, fi m = 0.27, and 
Ho — 71 krns" 1 Mpc -1 . All quoted magnitudes are on the 
Vega system. 



2 OBSERVATIONS, REDUCTION AND 
REDSHIFTS 

2.1 Sample Selection 

The objective of this study is to test the link between SMGs 
and QSOs at the era where these populations peaked. We 
have therefore selected radio-quiet QSOs with a redshift 
distri bution matched to SMGs ((z) ~2.3, IChapman et al.l 
120051 ) and with 850 /im fluxes ^ 5 mjy. The Ssso ^ k, 5 mjy 
flux limit selects comparatively bright QSOs and SMGs 
and is imposed so that the sources would be sufficiently 
bright to be detected in CO. We choose submm-detected 
QSOs from three main samples: we select five QSOs from 
a submm ph o tomet ry survey of optically-bright QSOs by 
lOmont et all (|2003T ) - HS1002+4400, J140955.5+562827, 
J154359.3+535903, HS1611+4719, and J164914.9+530316, 
we select two QSOs from the sub mm photomet r y sur - 
vey of X-ray absorbed QSOs by IStevens et alj (120051 ) 
- RXJ121803. 82+470854.6 and RXJ124913.86-055906.2, 
and finally we in clude two QSOs foun d in blank-field 
submm surveys bv lChapman et al.l ([20051 ) and G reve et al.l 
|2004l ). SMMJ131222. 35+423814.1 and MM J163655+4059, 
respectively, and SMM J123716. 01+620323. 3, an optically- 



selected QSO targett ed in submm photometry mode by 
Chapman et all (|2005h . 

The final sample comprises 10 sub mm-detected QSOs 
with similar fluxes to the SMG sample (|Greve et alj|2005l ; 
Z 5mJy at 8 50 ^m, assuming 8 50/xm/1.2mm colours typi- 
cal of SMGs; iGreve et aill2004l ). with a median of lOmJy. 
The SEDs of these QSOs demonstrate that the bulk of the 
mm/submm luminosity arises from dust emission and is very 
likely associated with star formation, r ather than an AGN 
jBeelen et al J 120061 ; lOmont et all 120031 ). Our sample spans 
six orders of magnitude in QSO i?-band brightness, enabling 
us to potentially trace how CO properties depend on QSO 
luminosity or Mbh (see Fig. [1} . 

2.2 Near-Infrared Spectroscopy 

The QSOs in this sample were spectroscopically identi- 
fied through rest-frame UV spectroscopy. However, since 
the rest-frame UV emission lines can be extremely broad 
(FWHM~ 5000 km s" 1 ), and can be offset in velocity with 
respect to the systemic redshift, precise systemic redshifts 
are required to ensure the CO emission falls within the PdBI 
receiver bandwidth. To obtain precise systemic redshifts for 
our targets we use the nebular emission lines of [Oiii]A5007 
and H/3 which are redshifted to near-infrared wavelengths. 

Seven of our submm-detected QSOs were observed 
with the United Kingdom InfraRe d Telescope (TJKIR T) 1- 
5 fim Imager Spectrometer (UIST; iHowatt et all 12003 ). We 
used the HK grism which provides a spectral resolution of 
A/AA ~ 1000 and covers a wavelength range 1.4-2.4/im, 
thus allowing us to cover H/3, [Oni] A4959,5007 and Ha si- 
multaneously for most of our QSOs. The QSOs were ob- 
served between 2006 Feb and 2006 Oct. All observations 
were taken in <; 0.8" seeing and clear conditions, and were 
carried out using the standard ABBA configuration where 
the QSO was nodded by 6-12" along the slit to achieve sky 
subtraction. Individual exposures were 240 s and the total 
integration times varied between 1.2 ks and 4.8 ks depend- 
ing on the Tf-band magnitude of the QSO (see Table [TJ. I n 
addition, three Q SOs in our sample h ave suitable archival 
spectroscopy from lTakata et all l|2006l) and lSwinbank et al.l 
l|2004l2006l ). 

The relevant orac-dr pipeline (|Cavanagh et al.ll2003l ) 
was used to sky-subtract, extract, wavelength calibrate, flat- 
field and flux calibrate the data (and in the case of integral 
field unit observations, forms the datacube). Spectra for the 
QSOs in our sample are shown in Fig. [2] To derive systemic 
redshifts, we simultaneously fit the H/3 and [Oni] A4959,5007 
emission lines with Gaussian profiles. The emission line flux 
ratio of the [Oni] A4959,5007 emission lines is fixed such that 
^4959/^5007 =2.9, and the [Oni] A4959,5007 line widths are as- 
sumed to be the same. The fit assumes a linear continuum 
plus Gaussian emission line profiles allowing the velocity 
centroid, flux and width of the H/3 line profile to be dif- 
ferent from [Oni] . This makes a total of seven parameters to 
be fitted. We also derive emission line fluxes and widths for 
the Ha emission lin e (which are import ant for deriving Mbh 
using the method of I Greene fc Ho 2005) in a similar manner, 
fitting a single Gaussian emission line profile superimposed 
on a linear continuum. The resulting redshifts and fluxes are 
given in Table [2] These redshifts were subsequently used to 
target the CO emission with PdBI (see i ]2.3l) . 
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Table 1. Summary of the new near-infrared spectroscopy and mm- wave PdBI CO obser vations of submm-detected 2 < z < 3 QSOs, 
including nine new CO observations and one published object from iBeelen et al The CO t eX p corresponds to the on-sourcc 

integration time with the equivalent of a 6-element array. 



Source 


Ncar-IR Observations 






CO Obscrvations a 












Datc(s) 




Dates 




Noise per channel^ 


Beam 


FWHM 


Detected? 






(ks) 




(hrs) 


(mjy beam - 1 ) 










HS1002+4400 


2006 Mar 08 b 


1.8 


2006 Sep 16,17,19,21,22 


7.2 


1.4 


4.6" 


X 


3.3" 


Y 


RXJ121803. 82+470854. 6 






2001 Dec 20, 22, 24 


9.0 


1.2 


6.7" 


X 


3.2" 


N 


SMM J123716. 01+620323. 3 


2005 Fob 16 d 


4.0 


2007 May 12,18,30, Jun 3,5,18 


16.0 


0.5 


4.7" 


X 


4.1" 


N 


RXJ124913. 86-055906. 2 


2006 Feb 25 b 


4.8 


2001 Aug 5, 2002 Apr 13,14,25, May 6 


9.4 


1.8 


8.8" 


X 


6.6" 


Y 


SMM J131222. 35+423814.1 


2004 Jun 25 d 


2.0 


2007 Jun 2,4,9,10,22,23 


9.3 


0.5 


5.2" 


X 


4.1" 


Y 


J140955. 5+562827 


2006 Mar 06 6 


1.2 


2002 May, June, 2003 Jan-Mar' 


37.0 


0.6 


2.4" 


X 


1.7" 


Y 


J154359. 3+535903 


2006 Feb 25 b 


4.8 


2006 Jun 4,5,9,11,19,26 


12.1 


1.2 


5.6" 


X 


3.5" 


Y 


HS1611+4719 


2006 Aug 17 b 


3.6 


2006 Jun 19,21,27, Jul 7 


6.3 


2.1 


5.0" 


X 


4.1" 


Y 


MM J 163655+4059 


2003 Aug 03 c , 2004 Apr 07 d 


2.4 


2007 Apr 30, May 21 


9.5 


0.4 


5.3" 


X 


4.4" 


N 


J164914. 9+530316 


2006 Mar 06 b 


1.8 


2006 Jul 20,23, Aug 15,18 


6.5 


1.6 


5.5" 


X 


3.7" 


N 



Notes: 

(a) All observed in the CO(3-2) transition, except for RXJ121803. 82+470854.6 observed in CO(2-l) 

(b) From UIST spectroscopy (this work) 

(c) From ISwinbank et all < 20041) 

(d) From Takata et al. (200G) and Alexander et al. (2008) 
(c) From .Bc clen et al ! M2004h 

(f ) With a channel width of 20 MHz 



Due to the very weak [Om]A5007 emission in 
RXJ124913.86, it was very difficult to measure a reliable 
redshift for the CO observations, and thus a range of possi- 
ble redshifts were covered on this QSO (z ~ 2.24-2.26). 

During the course of assembling the data for this 
project, we obtained a CO detection of HS1611 after se- 
curing a systemic redshift from identification of [Oni] with 
UIST. The positional centroid of the CO was offset by 
~1.3arcsec north of the optical and near-infrared QSO po- 
sition, possibly indicating a companion, and HS1611 was 
re-observed with the UIST integral field unit (IFU) for a 
total of 3.6 ks on 2006 Sep 01. The UIST IFU uses an im- 
age sheer to take a 3.3" x 6.0" field and divides it into 14 
slices of width 0.24". The dispersed spectra from the slices 
are reformatted on the detector to provide two-dimensional 
spectro-imaging, in our case also using the HK grism. We 
reduced the data using the relevant ORAC-DR pipeline which 
extracts, flatfields, wavelength calibrates the data and forms 
the datacube. We discuss these data in £|3. 1 .81 

2.3 Millimetre Interferometry 

We surveyed seven QSOs for CO emission using IRAM PdBI 
between 2006 June to 2007 June and two QSOs between 
2001 December and 2002 April in conditions with good at- 
mospheric phase stability and reasonable transparency. The 
observations were carried out using 3, 4, 5, and 6 anten- 
nae in the D configuration, giving a total of 3, 6, 10, and 
15 baselines, respectively. For uniformity with the previous 
observations observations of SMGs, which will comprise our 
comparison sample, we have aimed t o achieve the sam e 1- 
a sensitivity as for the SMC survey l|Greve et alj|2005h for 
the new CO observations of our sample of submm-detected 
QSOs: Sco An ~ 0.3 Jykms" 1 . We also includ e results on a 
single QSO from the literature, J140955.5 from lBeelen et ail 
(2004), to yield a final sample of 10 QSOs. The PdBI obser- 
vations are summarised in Table [T] 

For data taken prior to 2007 January, the spectral corre- 
lator was adjusted to detect the line with a frequency resolu- 



tion of 2.5 MHz at a bandwidth of 580 MHz. For data taken 
after January 2007, the spectral correlator was adjusted to 
detect the line with a frequency resolution of 2.5 MHz in the 
1 GHz band of the new generation receivers. The visibilities 
were resampled to a frequency resolution of 5 MHz. 

The data were calibrated, mapp ed and analyzed in 
the I RAM gildajO] software package l)Guilloteau fc Lucas! 
2000). Anomalous and high phase-noise visibilities were 
flagged. Passband calibrations were typically made using 
a nearby bright quasar. Phase and amplitude calibration 
within each observation were calibrated using frequent ob- 
servations of nearby quasars approximately every 20 min. 
The flux calibration was done using standard calibrators, in- 
cluding 3C273, MWC349, 1637+574, 1749+096, 0923+392, 
1243-072, and CRL618, for example. We estimate a conser- 
vative error on the flux calibration of 10-15%. 

For each QSO the central observing frequency of the 
3.1mm band receiver was tuned to either the redshifted 
CO(3-2) or CO(2-l) rotational transition, depending on 
the systemic near-infrared redshift derived from our spec- 
troscopy given in Table [2] In the case of RXJ124913.86, 
the uncertainty on the redshift was large enough that we 
performed two separate observations with different central 
frequency tunings of 106.160 and 106.720 GHz to ensure 
that the CO (3-2) line was covered. We have combined these 
data to cover more effective bandwidth and have centered 
the resulting spectrum at the mean central observed fre- 
quency, 106.435 GHz, corresponding to CO(3-2) redshifted 
to z ~ 2.2487. 

Naturally weighted datacubes were made and inspected 
for line emission close to the QSO position. Fig.[3]shows the 
velocity-integrated mm emission maps over the line in each 
field. In Fig. [4] we show the spectra of the CO emission in 
the brightest pixel of the central source in each velocity- 
integrated mm emission map, centered on the systemic red- 
shift of each system. We have included the maps and spec- 
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Figure 2. Restframe UV— near-infrared spectra for our sample of submm-detected QSOs. The insets show the region around the H/3 
and [Oni]A5007 emission lines (all insets cover the wavelength range 4675— 5100A), except in the case of J163655, and the redshifts 
derived from the [Oin]A5007 emission are indicated. The redshift for J163655 was taken from the redshifted Ha due to ambiguity in the 
[Oin]A5007 emission lying near the edge of the wavelength coverage of the spectrograph; the inset therefore shows the region around Ha. 
The region around 5100— 5500A in J1409 has been masked due to strong atmospheric absorption. 



tra of non-detections and in the latter cases we show the 
maps averaged over the central 500 km s -1 of the spectra 
(the average line width of our sample) and spectra extracted 
at the phase tracking centre (or at the continuum peak 
for J164914.9). In all cases, we ignore the small correction 
for the primary beam attenuation, given that the detected 
sources are all close to the phase centres of the maps. 



3 ANALYSIS AND RESULTS 
3.1 CO properties 

We present new CO data for nine QSOs, including five new 
CO detections. Four of our QSOs have CO detections with 
S/N^ 5, and an additional QSO is marginally detected in 
CO with S/N~ 2.5 with emission coincident with the phase 
centre (see Fig. [4j| . None of the QSOs appear to be resolved 
with the ~ 6arcsec beams (~ 50kpc FWHM). We include 



a previous CO detection reported bv lBeelen et al.l (|2004t ) in 
our sample of submm-detected QSOs. In total therefore this 
gives six detections in our sample of 10 QSO s, yielding a 
detect ion rate comparable to that for SMGs bv lGreve et al.l 
l|2005h of ~ 60%. 

Looking at the non-detections, it is unlikely that the 
cause of the failure to detect CO emission in these four QSOs 
is due to the systemic redshifts being wrong, as the quality of 
the optical/near- infrared spectral data for these is high and 
similar to the successful detections (see Fig. Similarly, it 
is unlikely that the velocity offset of the CO emission and the 
systemic redshift are larger than the correlator coverage (the 
rms deviation of znir and zco for our CO-detected QSOs is 
0.0046), since only a few syste ms detected at high-r edshift 
so far show larger offsets (e.g. lHainline et al.l [20041 . It is 
also unlikely that the CO line widths are 3> 1000 km s -1 
since our typical detected lines are < 550 kms -1 which is 
much less than the correlator bandwidth. Although extreme 



6 Coppin et al. 



CO line widths could be one reason for some of the non- 
detections, the most likely reason appears to be that the 
QSOs are simply too faint in CO to be detected given the 
depth of our observations. However, the non-detections do 
still provide useful upper limits on the CO luminosity and 
gas masses and therefore we include them in our analysis. 

Other than one QSO (J164914.9; £13.1.100 , we do not 
detect the continuum at 3.1 mm in any of the QSOs. This is 
consistent with the measured submm/mm fluxes and grey- 
body dust emission with a spectral form S„ oc v 2+l3 in the 
Rayleigh- Jeans regime of the spectral energy distribution 
with a dust emissivity /3 = 1.5. 

The luminosity, velocity width and spatial extent of the 
CO line emission can be used to place limits on the gas and 
dynamical mass of each system. For each CO-detected QSO 
a Gaussian profile is fit to either 20- or 30-MHz binned data 
by minimising the % 2 statistic. The 1 a errors on the best- 
fitting parameters are quoted where appropriate. To calcu- 
late the velocity-integrated flux density of the CO emission, 
we integrate the best-fitting Gaussian profile (note that we 
do not fit the line and continuum simultaneously) . The error 
on the velocity-integrated flux density is calculated by boot- 
strapping from a Normal distribution centered on the best 
fitting amplitude above, with a width given by the 1 a error 
in the fitted amplitude, while holding the FWHM and cen- 
tral velocity of the best-fitting Gaussian profile fixed. The 
noise in the averaged channel maps, constructed from the av- 
erage emission over the full width at zero intensity (FWZ1) 
in each QSO, is determined by summing over the weights 
and takes into account the phase noise visibilities. 

For the non-detections we calculate 3<r upper lim- 
its Sco = 3 a (Sv Avfwhm) 1 ^ 2 , where a is the channel- 
to-channel rms noise, Sv is the veloc ity resolution and 
Aitwhm is the line width (see e.g. iGreve et alj 120051 ; 
ISeaquist. Ivison fc Halll Il995r i. We use the channel-to- 
channel rms noise of the spectra binned to a resolution of 
100 MHz and adopt a line width of 500 km s -1 . Since our de- 
tections have typical line widths < 550km s -1 , this is a con- 
servative assumption. In addition, this line width assump- 
tion facilitates the dir ect comparison of o ur non-detections 
with those of SMGs in lGreve et~ai1 (120051 ). 

We have summarised the line properties of each source 
(the CO position, z C o, L' co(3 _ 2) and FWHM) in Tables [3] 
and|31 Offsets of the CO emission from the radio position are 
given in Table [3] a ssuming that the so urce is unresolved us- 
ing equation B2 in llvison et al.l (|2007t ) . and are less than 2 a 
significant in most cases. We note that we discover serendip- 
itous ^ 4<r sources in two of the fields surveyed which 
may be associated with gas-rich companions to the QSOs 
(SMMJ123716.01 and SMM J131222.35). O verdensities of 
submm sources have already been found by IStevens et al.l 
|2004l ) around similar QSOs on slightly bigger scales than 
those probed here. We now discuss the CO properties of 
each QSO detection and non-detection in detail, as well as 
any serendipitous detections. 

8.1.1 HS 1002+4400 

The integrated CO emission in this QSO is detected at 
~ 6 a, 0.8" nort h of the Sloan Digital Sky Survey (SDSS; 
iFan et all fl999h position (see Fig. |3). The CO(3-2) spec- 
trum of HS 1002+4400 has been binned to a frequency res- 



olution of 20 MHz (54 kms -1 ) and is shown in Fig. 0] The 
CO emission line is well fitted by a Gaussian profile with a 
FWHM of 640±160kms -1 , a velocity-integrated flux den- 
sity of 500(3-2) = 1-7 ± 0.3 Jykms -1 , and a CO redshift 
of zco = 2.1015 ± 0.0007, offset by 20 ± 65 kms -1 from 
the near-infrared systemic redshift. No significant contin- 
uum emission is detected from the line-free region down to 
a 1-a limit of 0.4 mjy. 



3.1.2 RXJ121803. 82+470854.6 

Neither CO (2-1) emission at the systemic redshift nor con- 
tinuum emission are seen in the map near the QSO posi- 
tion. To calculate a 3 a upper limit to the velocity-integrated 
CO(2-l) line flux of RXJ121803.82, we extract a spectrum 
at the QSO position (Table [2)| and rebin the spectrum to a 
frequency resolution of 100 MHz (360 kms -1 ), resulting in 
a channel rms of 0.5 mjy, and so derive an upper limit of 
0.6 Jy kms -1 , assuming a line width of 500 kms -1 , centred 
on the znir = 1.7416. No significant continuum emission is 
detected down to a 1-a limit of 0.2 mjy. 

We do find two ~ 4 a serendipitous detections in the 
datacube, offset in velocity from the central frequency tun- 
ing, and approximately 10-20 arcsec away from the phase 
centre. Neither show strong broad line detections in the spec- 
tra extracted at the peak locations, so we believe these are 
unlikely to be true sources. 



3.1.3 SMM J 12371 6. 01+620323.3 

Neither CO (3-2) emission nor continuum emission are seen 
in the datacube near the QSO position. To calculate a 3 a 
upper limit to the velocity-integrated CO(3-2) line flux, 
we extract a spectrum at the radio position (Table [2J) and 
rebin the spectrum to a frequency resolution of 100 MHz 
(270 kms -1 ), resulting in a channel rms of 0.2 mjy, giving 
an upper limit of 0.3 Jy kms -1 , assuming a line width of 
500 kms -1 at the systemic redshift from [Om]A5007 (Ta- 
ble [2j. No significant continuum emission is detected down 
to a 1-a limit of 0.2 mjy. 

We find a ~4<r detection in the datacube, offset in 
velocity from the central frequency tuning and approxi- 
mately 18 arcsec0 south of the nominal pointing position 
at 12 3715.36, 62 03 06.5 (J2000), indicating a redshift of 
zco = 2.0597 ± 0.0001, which is close to the QSO zsm, 
suggesting it may be a gas-rich companion of the QSO 
with a separation of 170 kpc. See Fig. [5] This position 
corres ponds to within .5 arcsec of an I — 24.7 counter- 
part JSmail et al.1 |2004|) with a radio flux of < 15 ^Jy 
IjBiggs fc Ivisonll2006r) which is also unde tected in the Chan- 
dra observations (|Alexander"e t al. 2003). The CO line has 
an integrated flux of 5*00(3-2) — 0.2 ± 0.06 Jy kms -1 and is 
very narrow (FWHM=100 ± 30 kms -1 ). 



2 Note that the map noise at this distance from the phase centre 
is only slightly degraded since the primary half power beam width 
(HPBW) at this frequency is ~ 45 arcsec. 
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Table 2. Summary of the near-infrared observed properties of the submm-detected QSOs. 
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Source 


Near 


-IR posit: 


ion (J2000) 


R 


K 


z 


NIB 


FWHMhq 


Sec 


FWHM [0 in] 






R. 


A. 


Dec. 


(Vc f 










(kmr 1 ) 


(xlO -19 W/m 2 ) 


(kms -1 ) 


HS1002+4400 


10 


05 


17.43 


43 46 09.3 


16.1 


14.2 


2.1015 


± 


0.0010 


10000 ± 1000 


750 ± 100 


1500 ± 300 


RXJ121803. 82+470854. 6 


12 


18 


04.54 


470851.0 


19.5 


17.0 


1. 


7416 








SMM J123716. 01+620323. 3 


12 


37 


16.00 


62 03 23.4 


20.2 


15.6 


2.0568 


± 


0.0013 


2100 ± 500 b 


48 ± 10 b 


1100 ± 200 


RXJ124913. 86-055906. 2 


12 


4!) 


13.85 


-05 59 19.4 


16.2 


13.5 


2.2400 


± 


0.0100 


4820 ± 200 


1000 ± 100 


700 ± 300 


SMM J131222. 35+423814.1 


13 


12 


22.32 


42 38 13.9 


20.2 


18.0 


2.5543 


± 


0.0010 


2600 ± 1000 b 


60 ± 10 b 


670 ± 300 


J140955. 5+562827 


14 


09 


55.50 


56 28 27.0 


17.0 


14.9 


2.5758 


± 


0.0050 


4250 ± 200 


220 ± 20 


600 ± 200 


J154359. 3+535903 


IS 


43 


59.44 


53 59 03.2 


16.8 


14.3 


2.3692 


± 


0.0015 


8280 ± 300 


770 ± 50 


1000 ± 200 


HS1011+4719 


16 


12 


39.90 


4711 57.0 


17.1 


15.1 


2.4030 


± 


0.0008 


4010 ± 300 


300 ± 20 


1000 ± 300 


MM J163655+4059 


16 


36 


55.79 


40 59 10.5 


23.2 


19.1 


2.6070 


± 


0.0006 


3000 ± 400 


16 ± 2 


1840 ± 400 


J164914. 9+530316 


16 


4!) 


14.90 


53 03 16.0 


16.9 


14.2 


2.2704 


± 


0.0009 


3100 ± 300 


580 ± 50 


1000 ± 200 



Notes: 

(a) All redshifts are derived from [Olll] , except for RXJ121803.82, where the rcdshift is derived from Civ, Cm] and Mgll by |Page et al.| j200li) using WHT 
data 

(b) Ha flux is calculated from the H/3 flux, assuming an intrinsic Shq/Sh^ ^3.1 \ Alexander et al 1120081) . and the given line width is for H/3 



Table 3. Summary of the mm and CO observed properties of the submm-detected QSOs. 



Source 


CO posit: 


ion (J2000) 


Offset 


Resolution 3 


z CO 


Sco A« h 


FWHMj!. 
(kmr 1 ) 


850 fim ffux 


1.2 mm flux 




RA 


Dec. 


(arcscc) 


(kpc) 




(Jykms- 1 ) 


(mJy) 


(mjy) 


HS1002+4400 


10 05 17.43 


43 46 10.1 


0.5 


33 


2.1015 ± 0.0007 


1.7 ± 0.3 


640 ± 160 




4.2 ± 0.8" 


RXJ121803. 82+470854. 6 








42 




< 0.6 




6.8 ± 1.2 b 




SMM J123716. 01+620323. 3 








37 




< 0.3 




5.3 ± 1.7 C 




RXJ124913. 86-055906. 2 


1249 13.91 


-05 59 20.1 


0.6 


125 


2.2470 ± 0.0016 


1.3 ± 0.4 


1090 ± 340 


7.2 ± 1.4 b 




SMM J131222. 35+423814.1 


13 12 22.23 


42 38 13.9 


0.9 


38 


2.5564 ± 0.0011 


0.4 ± 0.1 


550 ± 220 


3.0 ± 0.9° 




J140955.5+562827 e 


14 09 55.50 


56 28 26.4 




17 


2.5832 ± 0.0001 


2.3 ± 0.2 


310 ± 30 




10.7 ± 0.6" 


J154359. 3+535903 


15 43 59.43 


53 59 03.4 


0.7 


38 


2.3698 ± 0.0006 


1.0 ± 0.2 


520 ± 140 




3.8 ± 1.1° 


HS1611+4719 


16 12 39.90 


47 11 58.3 


0.6 


37 


2.3961 ± 0.0002 


1.7 ± 0.3 


230 ± 40 




4.6 ± 0.7" 


MM .7163655+4059 








39 




< 0.2 






2.2 ± 0.6 d 


J164914.9+530316 / 








39 




< 0.8 






4.6 ± 0.8" 



Notes: 

(aj 1.2 mm flux from lOmont et al I j2003D 

(b) 850 /im flux from jpage et al.| i200.l l) 

(c) 850 /im flux from Chapman et al. (20(35) 

(d) 1.2mm flux from lGreve et alj j2004) 

(e) CO parameters from ^eelenetal ■I J2OO4I) 

(f) Detected in continuum only 

(g) Approximate CO Resolution at the target redshifts, given the beamsizes in Table [T] 

(h) This quantity was obtained by fitting a Gaussian distribution to the CO spectrum (see text) 



3.I.4 RX J 12491 3. 86-055906.2 

The integrated CO emission is detected at ~ 5 a, l.larcsec 
southeast of the QSO position (Table [2]). The CO (3-2) spec- 
trum of RXJ124913.86 has been binned to a frequency res- 
olution of 30 MHz (84 kms -1 ) and is shown in Fig. Q] We 
have combined two sets of data with different central fre- 
quency tunings (see and have arbitrarily set the cen- 
tral frequency to 106.435 GHz, corresponding to CO(3-2) 
redshifted to z = 2.2487. The CO emission line is fitted 
by a Gaussian profile with a FWHM of 1090±340 kms" 1 , 
a velocity-integrated flux density of 500(3-2) = 1-3 ± 
0.4 Jykms- 1 , and a CO redshift of z CO = 2.2470 ± 0.0016, 
offset by -170 ± 140 kms -1 from the near-infrared redshift. 
This is our broadest CO line detection. No significant con- 
tinuum emission is detected from the line-free region down 
to a 1-a limit of 0.8 mjy. 



tion of 20 MHz (62 kms" 1 ) and is shown in Fig. H The 
CO emission line is fitted by a Gaussian with a FWHM 
of 550±220 kms -1 , a velocity-integrated flux density of 
•Sco(3-2) = 0.4±0.1Jykms _1 (S/N~3), and a CO redshift 
of zco = 2.5564 ± 0.0011 which is consistent with the near- 
infrared redshift within the 1 a errors. No significant contin- 
uum emission is detected from the line-free region down to 
a 1-a limit of 0.1 mjy. 

We find a ~ 4 a serendipitous emission line in the 
datacube at 1312 23.63, 42 3819.33 (J2000), approximately 
15arcsec northeast of the phase centre, corresponding to 
a low S/N detected line with a FWHM=450±150 kms" 1 , 
zco = 2.5408±0.0007, and S C o(3-2) = 0.40±0.17 Jy kms -1 
(S/N~ 2.5). See Fig. [5] We do not find an optical counter- 
part for this source down to K ~ 20 and z ~ 24. If real, this 
source is offset by 120 kpc and 1300 km s -1 from the QSO. 



3.1.5 SMM J '131222.35+423814.1 

The integrated CO emission of this QSO is marginally de- 
tected at ~ 2.5 a, coincident with the phase centre of the 
PdBI observations (Table |3}. The CO(3-2) spectrum of 
SMM J131222.35 has been binned to a frequency resolu- 



3 Note that the map noise at this distance from the phase centre 
is only slightly degraded since the primary half power beam width 
(HPBW) at this frequency is ~ 52 arcsec. 
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3.1.6 J 140955. 5+562827 

We briefly summarise the details for J140955.5 here since its 
detec ti on is published in llHainline et al 1 120041 : IBeelen" et al.l 
120041 ). IBeelen et all \200q ) report a CO(3-2) emission de- 
tection at the central observed frequency (corresponding 
to zco = 2.5832 ± 0.0001) that is well fitted by a Gaus- 
sian with a FWHM = 31 1 ±28 km s -1 , yielding S C 0(3-2) = 
2.3 ± 0.2 Jy kms -1 . 

3.1.7 J '1 54359. 3+535903 

We detect the integrated CO emission of this QSO at ~ 5 a, 
and coincides with the phase centre of the PdBI observations 
(Tabled Fig. g}. The CO(3-2) spectrum of J154349.3 has 
been binned to a frequency resolution of 20 MHz (58 km s _1 ) 
and is shown in Fig. [4] The CO emission line is well fitted 
by a Gaussian with a FWHM of 500±130kms~\ a velocity- 
integrated flux density of Sco(3-2) = 1.0±0.2 Jy kms -1 , and 
a CO redshift of z C o = 2.3698 ± 0.0006 consistent with the 
near-infrared redshift of the QSO within the 1 a errors. No 
significant continuum emission is detected from the line-free 
region down to a 1-a limit of 0.4 mjy. 

3.1.8 HS1611+4719 

The integrated CO emission in this QSO is detected at ~ 6 a, 
1.3 arcsec north of the phase center of the PdBI observations 
(see Fig. |3J. The CO(3-2) spectrum of HS1611+4719 has 
been binned to a frequency resolution of 20 MHz (59 km s~ 1 ) 
and is shown in Fig. 3] The CO emission line is well fitted 
by a Gaussian with a FWHM of 230±40kms"\ a velocity- 
integrated flux density of 600(3-2) = 1-7 ± 0.3 Jykms -1 
(S/N~6), and a CO redshift of z C o = 2.3961 ± 0.0002. No 
significant continuum emission is detected from the line-free 
region down to a 1-a limit of 0.5 mjy. 

The apparent 1.3 arcsec spatial and ~ — 600 kms -1 ve- 
locity offsets between the CO and near-infrared emission 
might be naturally explained if the CO is associated with a 
companion galaxy. To test this we search the UIST IFU dat- 
acube for a companion galaxy to the north, but are not able 
to identify any strong line emitting galaxies to a flux limit of 
~2xlO -17 Wm -2 (co rresponding to a n Ha star-formation 
rate of ~40M Q yr -1 ; lKennicutt[|l998t ). We note that the 
near-infrared spectrum shows an extended blue wing on the 
[Oiii] emission line which is coincident with the redshift 
of the CO emission (see i]2.2|l . Higher resolution mm spec- 
troscopy and/or deeper near- infrared spectral imaging may 
be required to determine whether the spatial and velocity 
offsets arise due to a companion galaxy, an outflow from the 
QSO, or a merger. This is the only example of a potentially 
significant offset between the CO and near-infrared position 
and/or redshifts. 

3.1.9 MMJ163655+4059 

Neither CO (3-2) emission nor continuum emission are seen 
in the map near the QSO position. To calculate a 3-<r upper 
limit to the velocity-integrated CO (3-2) line flux, we extract 
a spectrum at the radio position (Table [5} and rebin the 
spectrum to a frequency resolution of 100 MHz (310 kms -1 ) 
centred on the near-infrared redshift (Table [2} , resulting in 



a channel rms of 0.2 mjy, and we find an upper limit of 
0.2 Jykms -1 , assuming a line width of 500 kms -1 . No sig- 
nificant continuum emission is detected down to a 1-a limit 
of 0.2 mjy. 

3.1.10 J164914. 9+530316 

No CO (3-2) emission line is seen in this QSO, although con- 
tinuum emission is detected in the averaged channel map 
with a flux of 1.6 ± 0.28 mjy (S/N~6) at a position of 
16 4914.853, 53 0316.35 (J2000) (see Fig.©. 

We have determined that the continuum is dominated 
by synchrotron or free-free emission, with a radio spectral 
index, a', of 0.1, where S v oc u a , based on 1. 4 and 5 GHz 
fluxe s of 820±20 and 910±80 /x Jy, respectively l|Petric et all 
2006) which predicts a 3 mm flux of 1.3 mjy, in agreement 
with our detection. Based on the 3 mm flux, we estimate 
that approximately a third of the 1.2 mm flux (1.2 mjy of 
the Si. 2mm = 4.6 ± 0.8 mjy) is due to synchrotron/free- 
free emission rather than dust emission. We calculate a 3 a 
upper limit to the continuum-subtracted velocity-integrated 
CO(3-2) line flux of 0.8 Jykms -1 , by rebinning the spec- 
trum to a frequency resolution of 100 MHz (280 kms -1 ), re- 
sulting in a channel rms of 0.7 mjy, and assuming a line 
width of 500 kms -1 centred on znir = 2.2704. 

3.2 CO Luminosities and Gas Masses 

For each submm-detected QSO, we calculate the line lumi- 
nosity and estimate the total cold gas mass, M gaB , from the 
integr ated CO line flux following ISolomon fc Vanden Bout] 
(2005). We assume a line luminosity ratio of rz2 = »*2i = 
r = ^ / co(3-2)/^'co(i-o) = 1 (i- e - a constant bright- 
ness temperature). For our QSO sample, we derive a me- 
dian of L' C o(i-o) = (3-1 ± 0.9) x 10 10 Kkms -1 pc 2 for 
the non-detections and a median of £cc>(i_o) = ^~ 
1.2) x 10 10 Kkms -1 pc 2 for only the CO detections (see 
Table [3}- We then convert the CO luminosities into to- 
tal cold gas masses (including a correction factor for He- 
lium) using M gas = A/(H2+He) = a -^co(i-o) • We use 
q = O.8M0(Kkms -1 pc 2 ) -1 as the CO-to-gas conversion 
factor appropriate for local galaxy populations exhibiting 
similar levels of star formation activity to submm-bright 
galaxies or QSOs (e.g. ULIRG s). r — 1 and a — 0.8 are 
also used bv lGreve et all (|2005l ) for the SMG sample and so 
our QSOs can be directly compared with SMGs (if these val- 
ues are appropriate for both populations). The gas masses 
for our QSOs can be found in Table[4] The median gas mass 
of the entire sample is (2.5 ±0.7) x 10 10 Mq when the upper 
limits are included, and the median gas mass for the CO- 
detected QSOs is (3.4 ± 0.8) x 10 10 M . These gas masses 
are compar able to those observ ed in higher-redshift (z 4) 
QSOs (e.g. iRiechers et alj|2006h . 

Our observations are of J > 2 CO transitions, and it is 
possible that significant amounts of cold, possibly subther- 
mal molecular gas could be present, but only detectable in 
lower CO transitions. Moreover, if the gas is metal-poor, it 
is possible that the gas mass could be higher, or we could 
be missing clumpy dense g these CO observations pri- 
marily trace diffuse ISM. Fortu nately, these po s sibilit ies are 
unlikely to have a large effect: IRiechers et al] (|2006T > show 
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Figure 3. The velocity-integrated mm emission of our nine new submm-detected QSOs. The solid contours are 1, 2, 3... X ct (dashed 
lines represent equivalent negative contours), where <r refers to the noise in each velocity-integrated map (integrated over the width 
of the CO line) determined by summing over the weights and includes a contribution from the phase noise visibilities. These are the 
dirty maps (except for RXJ124913.86 which has been cleaned for presentation here to remove significant sidelobe structure) and the 
synthesized PdBI beams are shown in the insets as hatched ellipses for reference. Crosses indicate the radio or optical position of each 
source. Detections: The CO-detected sources shown here are HS1002, RXJ124913.86, J131222.35, J154359.3, and HS1611. Note that the 
CO detection of J131222.35 is marginal. J164914.9 shows a strong continuum detection but no CO detection (see ^3.1.101 . The positional 
uncertainties are ~ 0.5, 0.6, 0.9, 0.7, 0.6 and 0.5, respectively, for the CO or continuum detected sources (ranked in RA) based on the S/N 
of the detections and the beam sizes (see text). N on- detections: The three CO and continuum non-detections are shown for comparison 
and show the 3 mm emission over the central 500 km s^ 1 of the spectra (the average line width of our sample). 



that the CO emission in three well-studied z > 4 QSOs is 
well-described by a single centrally-concentrated molecular 
gas component and is highly excited. But a concern remains 
that if the CO-to-gas conversion factor for a faint extended 
component is higher (e.g. Galactic), a higher H2 mass may 
be hidden in such an extended component. Nevertheless, 
with suitable caution these observations can still be directly 
comparable to the SMG sample as they suffer from the same 
uncertainties and potential biases. 



3.3 Line Widths and Dynamical Masses 

The median CO line width of our sample is 550 ± 180 kms - 
(Table [3} . CO line widths can be directly converted into 
dynamical masses, assuming a size and inclination for the 
gas reservoirs. In the following we assume a disk model 
1 Solomon fc Vanden Boutl l2005t ). Following iTacconi et alj 
1 2006f ). we derive v^sin(i). where v c is the circular veloc- 
ity at the outer CO radius and i is the inclination of the 
gas disk relative to the sky plane, by dividing the FWHM 
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Figure 4. Millimetre spectra of the nine new submm-detected QSOs presented in our sample. The LSR velocity scale is relative to the 
near-infrared redshift (Table [2} of each QSO, except in the RXJ124913.86, where two central frequency tunings were used, reflecting the 
uncertainty of the systemic redshift of the source at the time of the CO observations were made, and so we have arbitrarily centered 
the spectrum on CO(3-2) redshifted to z = 2.2487. Detections: The CO-detectcd sources HS1002, J154359.3, J131222.35 and HS1611, 
have been binned to a frequency resolution of 20 MHz, and RXJ124913.86 has been binned to 30 MHz. The best-fitting CO redshift 
for each QSO is ob tained by fitting each spectrum with a single Gaussian distribution (overplotted) and is indicated by an arrow. See 
iBeelen et al. I d2004h for the CO spectrum of J140955.5. N on- detections: The CO-undetected spectra have been extracted from the phase 
tracking centre of the map and have been binned to a resolution of 20 MHz for comparison with the CO detections. The spectrum of 
J164914.9 has been continuum-subtracted, and the original spectrum is given as the dashed histogram. 



of the CO line by 2.4. We calculate the dynamical masses as 
A^dyn = Rv 2 esc 2 (i) / G ', where R is the radius of the gas disk 
which we assume to be 2 kpc, and list these in Table [4] With 
a typical resolution of ~ 40 kpc for our sample (see Table [3}, 
putting a 2 kpc radius limit is a major assumption (as are 
assumptions on the system inclination), although we explain 
why a 2 kpc disk is probably an appropriate assumption for 
our QSOs in We note that our estimates are conserva- 
tive as the dynamical masses will be h igher by a factor o f 
about two if a merger model is adopted (jGenzel et al.ll2003l) . 

For the CO-detected QSO sample we have a median 
dynamical mass of M(< 2 kpc) ~ (2.5 ± 1.6) x 10 10 csc 2 (i). 
The main uncertainties in the dynamical mass limits are the 
assumed disk size and the inclination angle, i. We note that 
one of our detected QSOs, RXJ124913, has a CO FWHM 
amongst the largest ever observed for a QSO at any red- 
shift, although our measurement is highly uncertain (the 
median FWHM for all CO-detected QSOs, including lensed 



and z > 4 objects, is 300 km s 1 ; ISolomon fc Vanden Bouti 
2005), suggesting that RXJ124913 is viewed at a high in- 
clination angle or that the CO disk size is larger than we 
have assumed. However, high resolution observations of the 
CO distribution are needed to better constrain the gas disk 
sizes, inclination angles, and hence dynamics of our QSOs. 

3.4 Far-infrared Luminosities, SFRs and SFEs 

We have derived Lfir, for the submm-detected QSOs from 
their restframe far-infrared fluxes and report these in Ta- 
ble [4] We scale a modified greybody model for far-infrared 
emission to match the 850 or 1200 /im photometry assum- 
ing a dust temperature of Td = 40 K and a dust emiss ivity 
factor of v 13 , with j3 = 1.5 (see e.g. ICoppin et alj|200sl ) and 
integrate the SED to obtain Lfir. Our assumptio ns for Td 
and f3 agree with the results of lBeelen et al.l (|2006l ) for three 
350 ^im-observed QSOs between redshifts of 1.8-2.6 with a 
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Figure 5. The velocity-integrated mm emission (left) and spectra (right) of serendipitous emission-line candidate 'companion' detections 
in two of the CO datacubes. See § I3.1,3l and l3.1,5l for details about the individual detections. Left: A ~ 4<r source appears in each map 
(solid contours are 1, 2, 3... X a), indicated by an asterisk, offset from the PdBI phase centre (indicated by a cross). See the caption of 
Fig. [3] for details. Right: 20-MHz binned spectra are shown for each serendipitous emission-line candidate companion source (see the 
caption of Fig. [4] for details). 



mean fitted Ta — (35 ±7) K when j3 is fixed to 1.6. Note that 
assum ing Td = 50 K effec tively increases Lfir by a factor of 
~2.5 ()Wang et al.ll20o"^ ). We have corrected the observed 
1.2 mm flux of J164914.9 by ~-1.5mJy, assuming a radio 
spectral index of 0.1 (see £13.1. 10[) for synchrotron emission, 
before computing Lfir. 

These Lfir es timates are comparable to those of SMGs 
i|Greve et al.ll2005l ) which is not surprising given the simi- 
lar submm flux selection criterion. These assume that Lfir 



is produced predominantly by st ar-formation as op posed 
to AGN activity (see discussion in lOmont et al, 2003). Re- 
cent observations have shown that the Lfir in QSOs is not 
generally contaminated by AGN, but is due t o star for- 
mation, even in the m ost powerful QSOs (e.g. iLutz et all 
120071; I Wang et alll2007l) . although there are clear exceptions 
(e.g. IWeiss et al.ll2007l ). 

This calculation yields a median of Lfir = (8.0 ± 1.9) x 
10 12 Lq for the submm-detected QSOs. We have converted 
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Lfir i nto a star-f o rmati on rate (SFR) for our sample fol- 
lowing iKennicuttl lll99Sl); SFR=1360 ± 320 M a yr" 1 . This 
relation assumes a Satpeter et al.l (|l955T j initial mass func- 
tion (IMF) and applies to starb ursts with ages less than 
100 Myr. See lOmont et al.1 (120011 ) for a discussion of the full 
range of conversion factors taking into account the uncer- 
tainty in the burst age, IMF, and metallicity. 

The star formation efficiency (SFE) is a measure of 
how effective a galaxy is at converting its gas into stars 
and can be represented by the continuum-to-line ratio, or 
Lfib./L' co , as this ratio presumably traces the star forma- 
tion rate per total amount of gas in a galaxy. We find a me- 
dian SFE=250±100L© (Kkms" 1 pc 2 )" 1 for the full sample 
of submm-detected QSOs. 



3.5 SMBH masses 



We de termine Mbh for our QSOs using the I Greene fc Hoi 
(|2005l ) virial Mbh estimator which calculates the BH mass 
from the Ha or H/3 emission line widths and fluxes. We can 
measure these lines from our near-infrared spectra (Fig. [2}: 
the observed emission line properties are given in Table [2] 
and the BH masses are given in Table [4] We find a me- 
dian Mbh of (1.8 ± 1.3) x 10 9 M Q . Note that six out of nine 
of the submm-detected QSOs have very large BH masses 
(Mbh Z 1O 9 M0) an d are among the larg e st masses known 
for local BHs (e.g. IValtonen et al.l 120071 ; iMacchetto et al.1 
1 19971 ; iHumphrev et al.l 120081 ). and that two objects in the 
sample are at the upper envelope of known BHs at any red- 
shift, confirming that these BHs are probably hosted by the 
progenitors of the most massive present-day ellipticals. 

Note that we have not corrected the Ha fluxes for 
extinction, since extinction corrections are rarely applied 
to broad-line fluxes when estimating virial BH masses, 
even when rest -f rame UV lines are used. Additionally, 
lAlexander et al.1 (|2008T ) do not find evidence for large 
amounts of extinction of the broad-line region in their sam- 
ple of broad-line SMGs (the mean extinction is only Av ~ 
1.2 mags, assuming a n intrinsic ratio of Ha/H/3=3.1 and a 
ICalzetti et al.l (|2000T ) reddening law). Applying this mean 
extinction correction to our sample would increase the BH 
masses by ~ 15% (0.1 dex). We also investigate whether the 
Eddington ratios, rj, for our submm-detected QSOs are sim- 
ilar to more typical opt ically selected QSOs in t he same 
redshift range. Following iMcLure fc Dunlod (I2004T ). we cal- 
culate the mass accretion rates and Eddington ratios, rj, for 
our submm-detected QSOs and find a range in rj consistent 
with typical QSOs in th e same redshift range (77 ~ 0.1-1; 
IMcLure fc Dunlod [2004) , with a median rj ~ 0.7, which is 
slightly higher than, but consistent with typical QSOs given 
the large uncertainties involved. 



4 COMPARISON OF THE 

SUBMM-DETECTED QSOS AND SMGS 

We now compare the CO properties of the submm-detected 
QSOs with SMGs. In order to fairly compare the SMG and 
QSO samples, we have discarded three strongly l ensed galax- 
ies fro m the CO-observed SMG sample from iGreve et al.l 
l|2005l ) that would not have been detected with ^ 5 mjy 



at 850 fim. These would lie below our sample flux limit ex- 
cept for the boost from the lensing magnification, and they 
thus probe an intrinsically fainter population. This leaves 
a sample of 17 SM Gs (of which 11 are CO-detected) from 
IGreve et al.l (|2005T ) with a median observed 850 /jm flux of 
8.2 mjy. 

We compare median values of the gas masses, CO 
linewidths, SFEs, Md yn , and Mbh/M sp i,. The error bars on 
these median values represent the spread of values, and have 
been obtained by taking the standard deviation of the me- 
dian values from 50 trials of randomly drawing out the ap- 
propriate number of values out of each population with re- 
placement. We note that the dominant uncertainties in M gas , 
Lfir, and Md yn are likely the assumed CO-to-gas conver- 
sion factor, the assumed j3 and Td in the SED fitting, and 
the assumed CO radius and inclination angle, respectively, 
and have not been included in each error budget, although 
the relative M gas , Lfir, and Md yn in our samples will be 
correct if the same assumptions hold for all our SMGs and 
submm-detected QSOs. 



4.1 The gas masses and star- formation efficiencies 
of submm-detected QSOs and SMGs 

The median gas mass of our CO-detected QSOs is M gas = 
(3.4 ± 0.8) x 10 10 M (see Fig. EJ). This is similar to the me- 
dian gas mass of the comparison sample of 11 CO-detected 
SMGs (M gas = (3.0 ± 0.5) x 10 10 M Q ). Including the up- 
per limits from the CO-undetected objects in each sam- 
ple yields median gas masses of (2.5 ± 0.7) x 10 10 and 
(2.4 ±0.5) x 10 10 M Q for submm-detected QSOs and SMGs, 
respectively. Not surprisingly, we find that the gas mass dis- 
tributions of both samples are statistically indistinguishable 
at the present sample sizes, as we show below. 

Given the moderate fraction of non-detections in 
both samples, we em ploy survival analysis (see e.g. 
Fcigclson fc Nelsonlll985l ) to estimate the intrinsic distribu- 
tions of the gas masses of submm-bright SMGs and QSOs. 
This includes the detections and upper limits, assuming that 
the data are censored at random (i.e. that the chance of 
only a gas upper limit being available for an object is inde- 
pendent of the true value of the gas mass). Given the cur- 
rent detection limit and small dynamic range in our sample, 
there is no strong apparent trend with the selection criteria 
(i.e. with submm flux) versus whether a particular source is 
detected or not, and thus the assumption of the data being 
randomly censored seems reasonab le. We employ the widely- 
used product-limit Av ni estimator (|Avni et al II 1980) ; sec also 
IWall fc Jenklr3l2004h to construct a maximum-likelihood- 
type reconstruction of the true distribution of gas masses 
in the QSO and SMG samples (see Fig. 0. We now test 
whether or not our censored samples of SMG and QSO 
gas masses are likely to have been drawn from the same 
dist ribution using the n on-pa rametric Gehan test, follow- 
ing |j^igelson^^Nelson| l|l985h . The Gehan statistic yields 
L n = 42 ±39 (equivalent to ~ Ict), revealing that the differ- 
ence between the gas mass distributions for SMGs and QSOs 
is about 35% likely due to chance. We consider 35% to be 
a high value, and therefore conclude that the gas masses of 
our QSO sample are indistinguishable from that of the SMG 
sample at the same epoch of z ~ 2-3 for the high mass end 
(where we are the most sensitive). 
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Table 4. Physical properties of our submm-detected QSO sample derived from the near-infrared and CO observations. 



Source 


L CO° 


T b 
FIR 


SFR C 


id, 


M dyn sin 2 


log M BH 




(xlO 10 Kkms- 1 pc 2 ) 


(X10 13 L Q ) 


(MQyr" 1 ) 


(XlO 1 ' Mq) 


(XlO 10 Mq) 


(M ) 


HS1002+4400 


4.2 ± 0.8 


1 1 + 13 
1 ' i -0.9 


+ 2200 
lauu -1600 


3.4 ± 0.7 


, ,+1.9 
a,a -1.8 


10.14 ± 0.2 


RX .1121803. 82+470854. 6 


< 2.4 


6 +0 ' 8 


1100+ 1400 


< 1.9 






SMM J123716. 01+620323. 3 


< 0.6 


5+ - 7 


„o n +1200 
aJU -620 


< 0.5 




8.1 ± 0.5 


RXJ124913. 86-055906. 2 


3.6 ± 1.0 


7 +0 ' 9 
"•'-0.6 




2.9 ± 0.8 


9 7+ 7 - 
a ' ' -5.1 


9.76 ± 0.2 


SMM J131222. 35+423814.1 


1.2 ± 0.4 


3+ ' 4 
U.J_ 2 


520+ 670 


1.0 ± 0.3 


2 5+ 2 ' 4 


8.2 ± 0.5 


VV96 J140955. 5+562827 


8.2 ± 0.6 


'-2.6 


4600+ 5000 


6.6 ± 0.5 


8+ ' 2 


9.28 ± 0.2 


VV96 J154359. 3+535903 


3.1 ± 0.7 


1 0+ 1 ' 3 
-0.7 


170Q + 2200 
"""-1200 


2.5 ± 0.6 


2 0+ 1 ' 2 


10.13 ± 0.15 


HS1611+4719 


5.1 ± 0.8 


1 2+ 1 ' 4 


2100 + 2400 


4.1 ± 0.6 


4+ ' 2 
u -*— 0.1 


9.26 ± 0.2 


MM J163655+4059 


< 0.8 


6 +0 ' 7 


„ 50 +1200 


< 0.6 




8.4 ± 0.5 


VV96 J164914. 9+530316 


< 2.2 


°' 8 -0.7 


1400+ 1600 


< 1.8 




9.16 ± 0.15 



Notes: 

(a) This quantity was obtained by fitting a Gaussian distribution to the CO spectrum (see text) 

(b) Derived assuming (3 — 1.5 and T — 40 K (sec text) 

(c) Derived following Kcnnicutt (1998) (sec text) 

(d) Derived assuming a — 0.8 

(c) Derived assuming a disk radius of 2 kpc 

(f) 1.2 mm flux is corrected for synchrotron emission before computing 
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tically significant difference in the SFE between the CO- 
detected and undetected QSOs in our small sample given the 
current detection limits. What is the cause then of a given 
submm-detected QSO in our sample, which all have similar 
Lfir, being either CO-detected versus CO-undetected? A 
possible explanation could be that the CO-detected submm- 
detected QSOs are physically larger and contain proportion- 
ally more massive CO reservoirs than the CO-undetected 
QSOs, although time-dependent star formation, differences 
in the SEDs or filling factors could also be causes for the 
scatter in the Lfir-Lcq correlation (see Fig. [7)l . 
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men 
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sample of SMGs extending this trend 
out to the highest far-infrared luminosities (> 1O 13 L0). For 
comparison, in Fig. [7]we have plotted our QSOs on the L^o - 
Lfir diagram along with LIRGs, ULIRGs and SMGs. The 
CO-observed SMGs and QSOs seem to lie along the relation 
within the considerable uncertainties in their far-infrared lu- 
minosities, suggesting that the SMGs and QSOs are under- 
going similar star formation modes, despite the QSOs host- 
ing luminous AGN activity. 



Figure 6. Distribution of the gas masses of the submm-bright 
QSO and SMG populations. We show an estimate of the cold 
gas mass distribution of submm-bright QSOs and SMGs, in- 
cluding the CO detections and the gas mass upper limits from 
the non-detections, constructed using the Avni estimator for 
censored data (see text). Note that the bins containing points 
logMgas < 10.2 and 10.4 for QSOs and SMGs, respectively, have 
been made arbitrarily large here since the Avni estimator requires 
> 1 detection in the lowest bin. The Gehan statistic reveals that 
the distributions are not significantly different. 



Within our QSO sample we find a broad similarity in 
the star formation efficiency of the CO-detected and unde- 
tected QSOs (see § I3.4|> . assuming conservatively that the 
CO-undetected subset lies just below the detection thresh- 
old. A Kolmogorov-Smirnov (KS) test reveals a 25% chance 
that the two subsamples are drawn from the same parent 
population, indicating that we have not detected a statis- 



To quantify this we compare the SFE of the QSOs to 
the SMGs, where Lfir for the SMGs has been calculated 
using the method outlin ed in i)3.4l and their L' co are taken 
from iGreve et all <|2005h . As Fig. [7] shows, the SMGs ap- 
pear to have similar SFE with respect to QSOs on average, 
although this probably reflects the similar submm flux se- 
lection criteria. Given that the median Lfir and M gas of the 
QSO and SMG samples are similar, (8.0±1.9) x 10 12 L and 
(7.9 ± 0.8) x 10 12 L Q and (2.5 ± 0.7) x 10 10 M Q and (2.4 ± 
0.5) x 10 10 M (including all CO-detected and -undetected 
QSOs and SMGs), respectively, the SFEs are also similar: 
medians of 250 ± 100 and 260 ± 50 L (Kkms -1 pc 2 ) -1 , for 
QSOs and SMGs respectively. We perform a Gehan test on 
these censored data sets and find L n = 56 ±35 (~ 1.6 a), re- 
vealing that the difference between the distribution of SFE 
in SMGs and QSOs is about 15% likely due to chance, indi- 
cating that we have not detected a significant difference in 
the SFE between our small samples of SMGs and QSOs. 
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Figure 7. Comparison of the CO and far-infrared luminosities for 
QSOs, LIRGs, ULIRGs, and SMGs. The line is the best-fitting 



ULIRGs and SMGs from lGreve et al.l 120051) . The CO-observed 
SMGs and QSOs appear to occupy the same part of the diagram 
(filled circles/stars represent CO-detections and open circles/stars 
are CO-nondetections), and lie just above the relation though the 
offset is not significant given the large errors in Z/fir- This is a 
particularly useful diagnostic, since it does not depend on the 
CO-to-gas conversion factor. 



4.2 The CO line profiles and dynamical masses of 
submm-detected QSOs and SMGs 

A comparison between the observed CO line profiles of 
SMGs and QSOs requires care as we know that there is an 
inherent difference in the types of lines seen in SMGs versus 
QSOs: five of the 11 SMGs detected in CO in ou r comparison 
sample show hints of do uble-peaked profiles ([Greve et al.l 
l2005l ; lTacconi et alJlioOrj ), while none of our QSO CO emis- 
sion lines appear double-peaked (which is perhaps not sur- 
prising given our small sample size and low SNR data). 
These double-peaked profiles could indicate merger activ- 
ity (as likely for SMGs from h igh-resolution obse rvations) , 
or a rotating gas disk/ring (e.g. lGenzel et aDfeOOSt ). In both 
cases, measuring the separation between the peaks is a more 
appropriate measure of the dynamics of the merger or ro- 
tating disk system than the FWHM of the individual peaks. 

In cases where double-peaked profiles are present we 
need to fit two Gaussian profiles simultaneously in order 
to derive the velocity offsets of the peaks. We have there- 
fore refit all of the SMGs and QSOs with single and double 
Gaussian profiles using the \ 2 statistic as a measure of the 
goodness-of-fit. We find that 4/11 of the CO-detected SMGs 
are better fit at the > 3.5 a level wi th a double Gaussi an pro- 
file over a single Gaussian profile (|Greve et al.ll2005h . None 
of the submm-detected QSOs are significantly better fit by 
double Gaussian profiles, although we note that the random 
chance of drawing six single peaked profiles (with replace- 



ment) out of a sample containing seven single-peaked and 
four double-peaked profiles (i.e. the SMG sample above) is 7 
percent. The fact that the QSO profiles appear to be single- 
peaked could indicate that our data are too low SNR to de- 
tect double-peaked profiles (especially for a merger, where 
both the width and peak strength of the lines can be differ- 
ent), or that we are viewing the QSOs almost face-orQ if the 
CO is contained in a disk, or that they are at a late stage of 
a merger, or that they are not in the course of a merger. 

For our six submm-detected QSOs we compare the CO 
FWHMs and dynamical masses with the CO-d etected sub- 
samp l e of 11 SMGs taken from the literature (|Greve et alj 
120051 ; iTacconi et all [20061 ') . The median FWHM of the six 
CO-detected QSOs and seven single-Gaussian profile fit 
SMGs are (550 ± 180) and (600 ± 130) kms -1 , respectively. 
When we include the double Gaussian profile fits in the 
median estimate, taking the separation between the peaks 
as a proxy for the dynamics of the system, the median 
for the SMGs slightly decreases to (530 ± 110) km s" 1 . We 
have plotted the histogram of the submm-detected QSO 
a nd SMG line widths f or comparison in Fig. [8] (c.f. fig. 2 
in lCarilli fc Wandl2006l ). 

The CO FWHM of QSOs and SMGs appear to be sim- 
ilar. To quantify this we use a standard two-sided KS test 
which reveals a 95% probability that the two distributions 
are drawn from the same parent population (this drops to 
77% if only the single-line Gau ssian profile SMGs are used). 
This is in stark contrast with ICarilli fc Wa ne (2006), who 
find only a 0.9% probability for an inhomogeneous sam- 
ple of 12 SMGs and 15 QSOs from a literature compi- 
lation (|Solomon fc Vanden Bout] [20051. However, we note 
that they have used the SMG single-Gaussian profile fit line 
widths even in cases where clearly double-peaked CO pro- 
files are apparent, leading to artificially broad FWHM in 
these cases. Their analysis is also based on a mixed sam- 
ple over a wide redshift range and includes strongly lensed 
sources. Also, the QSOs studied here have larger linewidths 
than have been reported for other QSOs and they are also 
those with the lowest SNRs and thus the most difficult to fit. 
Therefo re the difference s in ou r results for z ~ 2 QSOs and 
those in lCarilli fc Wand (|20061 ) for higher redshift QSOs can 
probably mainly be attributed to the small sample sizes in- 
volved and the somewhat different approaches (rather than a 
real statistically significant effect) when all possible sources 
of error are taken into account. 

Since the CO line widths of SMGs and QSOs are simi- 



4 An intrinsically double peaked profile could appear single- 
peaked if the velocity offset between the peaks is small enough 
(e.g. at high SNR a close merger could show an asymmetric line 
profile, hinting at two components). We attempt to put a limit 
on the inclination angles in our data by investigating a simple 
test case. We take two Gaussian distributions with a = 150— 
220kms~ 1 (similar to those found for double-peaked SMGs) and 
equal line fluxes and shift them to a velocity separation such 
that the FWHM of a single-Gaussian profile fit is ~ 550kms — 1 
(our QSO sample median FWHM), yielding typical velocity sep- 
arations of 300-150 km s . This translates into a conservative 
limit on the inclination of the system of i < 40° by using the 
relation v ^, s = v c - lrc sin(i) and taking v ^, s ~ 250kms — 1 and 
v c - 1IC = 400kms — 1 (the latter value is what is observed typically 
for SMGs by fitting double Gaussian profiles). 
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Figure 8. Comparison of the distributions of CO FWHM for 
11 SMGs and 6 QSO with detected CO emission. The distribu- 
tions have medians of (530 ± 110) km s — 1 (for a mix of single and 
double Gaussian profile fits), and (550 ± 180) kms" 1 for SMGs 
and QSOs, respectively. A KS test reveals a 95% probability that 
the SMG and QSO CO FWHM distributions are drawn from the 
same parent population (see text). The histograms have been off- 
set slightly for clarity. 



lar, it may follow that the dynamical masses of SMGs and 
QSOs are also similar. The uncertainties in this are if the gas 
is distributed similarly in both populations (i.e. either a disk 
or merger scenario), if their sizes (R) are the same, and if 
their average inclination angles (i) are similar. R ~ 2kpc 
is believed to be appropriate for the gas d i stribution i n 
SMGs (e.g. iGreve et all 120051 ; iTacconi et ail 120061 , l2008t) . 
O ne of our QSOs, J14 0955.5, has been observed at PdBI 
by iBeelen etafl (|2004T ) in CO(7-6) with a 1.0" x 0.5" res- 
olution and was unresolved, giving an upper limit to the 
gas reservoir of 5 kpc (though this could be too high a J 
transition to trace the bulk of the gas reservoir). No sizes 
have been measured for our sample of QSOs, although the 
CO emission seems to be confi ned to a compact re gion in 
other samples of high-z QSO (|Riechers et al.l l2006t ) as in- 
ferred by compar ing a range of CO line ratios (see also 
IWalter et al1l2004 who observe g as distributed over a 2.5 kpc 
radius and IRiechers et all 20081) . This agrees with the re- 
cent results from Maiolino et al.l (|2007l ) who find that CO 
emission in a z ~ 5 QSO is unresolved with a beam size 
of ~ 1 arcsec, implying that the molecular gas is contained 
with in a compact region w ith a radius of < 3 kpc. The QSO 
from lMaiolino et al.l (|2007r ) was observed in a high J transi- 
tion, CO(5-4), and one might worry that such high transi- 
tions might not trace the bulk of the reservoi r. CO (1-0) and 
(2-1) observations of other z z 4 QSOs (e.g. IRiechers et al.l 
2006) have confirmed that the high J transitions do trace the 
bulk of the reservoir, and thus provide good size constraints 
(this is in contrast to what is observed in most nearby galax- 



ies). Spatially resolved ob servations exist for l ocal QSOs, 
albeit the data are sparse: iDownes fc Solomon! J 19981) mea - 
sure a R > lkpc disk in Mrk 231; IStaguhn et alj (I2004T ) 
find a R > 1 kpc ring-like molecular gas distribution; and 
iKrips et al.l l|2007l ) estimate R ~ 3 ± 1 kpc for the CO emis- 
sion in a low luminosity QSO. Therefore, assuming 2-3 kpc 
radius regions will probably still be reasonable once these 
galaxies are resolved. 

While the average orientation of the SMGs is likely 
to be close to the value expected for random orientations, 
i ~ 30° , s ince they are not pron e to orientation selection ef- 
fects fsee lCarilli fc W ane 2006]), the submm-detected QSOs 
could be orientated preferentially with respect to the sky 
plane due to their selection as optically luminous QSOs in 
the first place. We find evidence for this from a comparison 
of the [Om] line width (see Table [2} , a proxy for the mea- 
sure of the stellar velocity dispersion (supporting the idea 
that the narrow-line region gas is in orb ital motion in the 
gravi t ational potential well of the bulge ; iNelson fc Whittle! 
1 19961 ; lNelsonll2000l ; iBonning et al.ll2005h . with the CO line 
width. We find that the [Om] FWHM are a median fac- 
tor of 1.9 ± 0.5 broader than the CO line profiles (with a 
large scatter; ratios range from 0.6-4.3), suggesting indeed 
that inclination effects co uld be important in the CO line 
widths of QSOs (see also IShields et al.ll2006T ). The lack of 
the double-peaked CO profiles in the submm-detected QSOs 
could be interpreted as a later phase of a merger in the typ- 
ical proposed evolutionary sequence, with SMGs represent- 
ing an earlier phase of the merger. We also know that at 
least two of our QSOs must have % < 20 (alternatively the 
assumed standard CO-to-H"2 conversion factor could be too 
high or else the assumed 2 kpc radius is too small) or else the 
cold gas mass estimates exceed the dynamical masses (even 
given the large error bars; see Table [4| . Conversely, some 
of our broadest CO line detections suggest that their incli- 
nation angles might be larger than 20°. Other constraints 
come from the obscured:unob scured AGN rat i o (~ 5-10) 
for a sample of SMGs, where lAlexander et al.l (|200Sf ) sug- 
gest that QSOs could be seen close to the line of sight 
(i = 18-25°). We thus adopt the best available estimates 
of the average inclination angles for SMGs and QSOs of 30° 
and 20°, respectively, in calculating the dynamical masses 
below. High-resolution CO observations could help to place 
tighter constraints on the inclination angles of our systems. 

Assuming a 2 kpc radius and an average inclination 
angle of 20° for the submm-detected QSOs and 30° for 
the SMGs, we find median dynamical masses, Md yn oc 
FWHMco 2 R esc 2 (i), of M(< 2 kpc) ~ (2.1 ± 1.4) x 10 11 
and (0.9 ± 0.4) x 10 11 Mq, respectively. This results in gas- 
to-dynamical mass fractions of ~ 15% and 30% for QSO and 
SMG populations, respectively. The gas mass accounts for 
a significant fraction of the dynamical mass of both classes, 
indicating that the host galaxies are at an early evolutionary 
stage. 

4.3 The BH-spheroid ratios of submm-detected 
QSOs 

We place a constraint on the spheroid mass by taking the 
dynamical mass estimate from H4.2I assuming i = 20° and 
a spheroid radius of 2 kpc or less (as we assumed for the 
CO emission), yielding a median spheroid mass of M 3p h = 
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Figure 9. The variation in Mbh/A^spIi with redshift, z, for differ- 
ent populations of galaxies and AGN. W e show the ratio meas ured 
for BHs in local spheroids (dotted line; lHarine fc Rix 20 04|) and 
those inferred for high-redshift radio-loud AGN I McLure_et ii af 
20061) . QSOs l|Peng et al.l I2006T ). and SMGs ^Alexander et al. 
2008 ). We compare these to the ratio determined for our CO- 
and submm-detectcd QSOs . We find goo d agreement between our 
estimate and that of Pen g et al. (2006) for similar, optically lu- 
minous QSOs at z ~ 2, supporting the assumptions used in our 
analysis, although we show how our estimate changes for different 
assumed inclination angles (dashed lines). Similarly, the typical 
QSO detected in our survey has a significantly higher A/bh/ M S ph 
ratio than inferred for z ~ 2 SMGs. 



(2.1 ± 1.4) x 10 11 M . In principle one should subtract the 
gas and BH masses, although we note that doing so makes 
a negligible difference in general, given the large error bars 
on Mdyn, and we also have an additional complication that 
two of our CO-detected QSOs have Af gas > AI dyn . Other 
uncertainties to our M ap h estimates include the possibility 
that the cold gas in submm-detected QSOs is more widely 
distributed (e.g. due to entrainment in outflows) than we 
have assumed here (^ 2kpc), which would mean that all of 
our dynamical (and hence spheroidal) masses have been un- 
derestimated. Similarly, if the spheroid is significantly more 
extended than the CO emission then the spheroid masses 
will be larger. For example, assuming a characteristic radius 
of ~ 10 kpc for the spheroid would result in an increase in 
the mass of a factor of ~ 5 times. 

Based on these assumptions we estimate MBn/M sp h ~ 
(Slg 1 ) x 10 -3 for our QSO sample. This is approximately an 
order of magnitude larger than the local ra tio of (1.4±0.4) x 
10~ 3 for the BHs in local s pheroids (Fig. l9l [Marconi fc Hunt] 
120031 ; lHaring fc Rixll2004r i. Thus our CO survey of z ~ 2 
QSOs appears to support the claim that M B H/M BD h for 
z > 1 AGN is higher than observed locally (e.g. IPeng et al.l 
2006; iMcLure et all [20061 ). Our M BH /M sph ratio could be 
made to match the local relation by either decreasing the 
inclination angles to ~ 6° or by adopting a spheroid extent 



ten tim es larger than we have assumed. Our results and 
those of IPeng et alj (|2006l ) suggest that BH growth may oc- 
cur more rapidly than spheroid formation in QSOs at z ~ 2 
(although both studies could be prone to strong selection ef- 
fects; see I Alexander et al.ll2008T ). Similar CO-based studies 
of rarer individ ual QSOs at higher redshifts also support this 
conclusion (e.g. lWalter et al.ll20o3 ; iMaiolino et al.ll2007l ). 

This result contrasts strongly with the MBH/M sp h 
claimed for SMGs, whi ch lie just below the local ratio (Fig. [^1 
I Alexander et al]|2008l ). However, as mentioned in H3.5I six 
out of nine of our submm-detected QSOs have very mas- 
sive BHs (A/bh Z 10 9 Mq ) which are far more massive 
than those inferred for typical SMGs ([Alexander et al]|2008l ; 
M B h ~ 10 8 Mq ) and would place them among the largest 
known masses for local BHs. They are therefore likely to 
be considerably rarer than typical SMGs. To make a fairer 
comparison between SMGs and submm-detected QSOs we 
need to focus on those CO-detected submm-detected QSOs 
with M B h more similar to that found for SMGs. We have 
three submm-detected QSOs in our sample with BHs more 
similar to those of SMGs, unfortunately two of these are 
undetected in CO, so we cannot determine their dynamical 
masses. Nevertheless, the one example with a marginal CO 
detection, SMM J131222.35, has an implied spheroid mass 
of M sph = 2.1 x 10 11 Mq (assuming i = 20°) and hence 
Mbh/Mspi, ~ 9 x 10" 4 . This is an order of magnitude be- 
low the ratio for the more optically luminous QSOs and 
places SMM J131222.35 approximately on the local relation 
(Fig. [9}. While based on a single source marginally detected 
in CO, this result hints that submm-detected QSOs with 
Mbh ^ 10* M might be the SMG-QSO 'transition' ob- 
jects we are interested in. We discuss this possibility further 
in 13 



5 DISCUSSION OF THE EVOLUTIONARY 
STATUS OF SUBMM-DETECTED QSOS 

In this section we discuss the various observations which 
constrain the evolutionary relationship between SMGs and 
submm-detected QSOs. For the purposes of this discus- 
sion, we ada p t the proposed evolutionary pathway of 
ISanders et al.l l|l988l ) and so test the scenario where the 
SMG population evolves through a submm-detected QSO 
phase, into a submm-undetected QSO and finally a passive 
elliptical. Now we ask - are the observed properties of the 
SMGs and submm-detected QSOs consistent with their ob- 
served gas depletion, amount of BH growth, and space den- 
sities required to link them in this evolutionary sequence? 

To provide a characteristic timescale for our discussion 
we first determine the lifetime of the submm-detected QSO 
phase. We assume that the far-infrared luminosities of the 
submm-detected QSOs arise from dust-obscured star forma- 
tion and that the molecular gas reservoirs we detect via their 
CO emission are the fuel for this star formation. In this way 
we determine that the QSOs have enough cold gas to sus- 
tain their current episode of star formation (assuming a star 
formation efficiency of 100%) for repletion ~ M gas /SFR ~ 
2.5 x lO lo M /136OM yr- 1 ~ 20Myr. Tdepletion is On the 
order of a couple of dynamical times within our assumed 
CO radius of 2 kpc (rd yn = R/v c < 9 Myr), which could be 
defined as a 'maximal starburst' (e.g. iTacconi et al.ll2006l ) 
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as is commonly found with other populations of starburst 
galaxies such as ULIRGs at low-redshift and SMGs and 
other QSOs at high-z. Since we are most likely catch- 
ing the QSOs halfway through the current star formation 
episode, we infer a submm-detected QSO phase lifetime of 
2 x r dc piction = 40 Myr. The reader should of course note that 
Tdepietion and our estimated lifetime are both lower limits if 
Lfir includes a significant contribution from AGN activity, 
or if the starburst is discontinuous, or if the IMF is more top 
heavy than Salpeter's IMF, or if the gas reservoirs con tain 
significant masses of sub-thermal gas dGreve et alj|2006h . or 
if the CO reservoir is widely distributed on scales S> 30 kpc 
(not that there is any evidence that such galaxies exist). 

Next we note that in our submm-detected QSO sam- 
ple CO is typically detected in only the optically luminous 
QSOs (see Fig. [I}, i.e. those QSOs with the most massive 
BHs (A/ BH > 10 9 M Q ; see Fig. OH). As we discussed in P31 
these BHs are much more massive than found in typical 
SMGs - is it feasible for the BHs to grow sufficiently be- 
tween these two phases? In Fig. [TU]we plot Af gas as a func- 
tion of Mbh for our sample of s ubmm-detected QSOs and 
for an 'aver age' SMG (Mbh fr om I Alexander et al.ll200Sl and 
M gas from iGreve et all l2005h . We indicate gas consump- 
tion and BH growth evolutionary tracks for the average 
SMG on timescales of 10 and 20 Myr, assuming a SFR of 
lOOOM yr _1 and Eddington-limited BH growth. Fig. \TU\ 
demonstrates that an SMG could evolve into a submm- 
detected QSO with M BH ^ 10 8 M Q and A'/ gas ~ 0.6 x 
10 10 M Q (the median M gas of the CO-undetected QSOs) 
in a reasonable timescale, whereas an average SMG would 
need substantially (~ 10 times) more gas and more time (;a 
SMG lifetime) to evolve into the optically luminous submm- 
detected QSOs in our sample with Mbh — 4 x 10 9 Mq and 
M gas ^ 3.0 x 1O 1O M . 

The gas consumption and BH growth timescales sug- 
gest that the optically luminous submm-detected QSOs in 
our sample cannot be the direct descendents of typical 
SMGs. Indeed the CO-detected optically luminous QSOs 
in our sample are extremely luminous AGN (Mb ~ —28, 
lOmont et all 120031 ) with Mbh > 10 9 M Q and will most 
likely evolve into the most massive elliptical galaxies found 
in local rich clust ers. These QSOs are so rare (~ 2 deg -2 
over all redshifts; iFan et all l200ll) that they are unlikely 
to be associated with any known SMGs, since Sssopm > 
4.5 mjy SMGs at z = 1-3 have much higher surface den- 
sities (~ 600de g~ 2 ; ICoppin et al.l I2006T) or space densities 
(~ 10" 5 Mpc- 3 ; IChapman et alj|2005l ). Even correcting for 
the relative timescales of the two phases still leaves an or- 
der of magnitude difference in the space densities. Based on 
the relative number or space densities, these optically lumi- 
nous submm-detected QSOs are unlikely to be involved in 
the SMG/QSO evolutionary sequence we are interested in 
testing. 

In contrast, the less lumino us submm-det e cted QSOs 
appear to be consistent with the ISanders et al.l (|l988h pro- 
posed evolutionary sequence, linking them to typical SMGs 
(see Fig. 1 1 Op . Thus we concentrate on the three submm- 
detected QSOs (including one marginal CO detection and 
two upper limits: SMM J131222.35, SMM J123716.01 and 
MMJ163655) with BHs more similar to that of SMGs of 
Mbh — 10 8 M . The median gas mass and SFR of this 
subsample are 0.6 x 10 10 Mq and 850 Mq yr _1 , respectively. 



Hence, these systems can sustain their current star forma- 
tion episode for Tdepietion ;£ 7 Myr or a total lifetime of 
the submm-detected QSO phase of g 15 Myrs. This esti- 
mate of the lifetime of submm-detected QSOs can then be 
used to relate them to SMGs as follows. The incidence of 
blank-field submm so urces identified with QSOs is ~ 4% 
|Chapman et alj|2005l ). In our test scenario this is the frac- 
tion of an SMG lifetime in which it appears as an submm- 
detected QSO. Conversely, about 15-25% of optically lu- 
mino us QSOs are detected in the submm l|Omont et alj 
l2003h . Assuming that all bright SMGs go through a QSO 
phase, and that all QSOs are formed via SMGs and there 
is no strong mass or luminosity dependence of their life- 
times, then these relative fractions imply that the total 
QSO lifetime (i.e. submm-detected plus submm-undetected) 
is about 20%/4% = 5 times shorter than the average lifetime 
of an SMG. Taking our estimated lifetime of less luminous 
submm-detected QSOs, <; 15 Myrs, these detection fractions 
imply total lifetimes for QSOs of % 70 Myrs and for SMGs 
g 350 Myrs. These estimates are consistent with indepen- 
dent estimates of the lifetimes for QSOs of ~ 20 -40 Myr 
|Martini fc Weinberg! lioOll ; iGoncalves et all 120081) and ~ 
100-3 00 Myr for SMGs from modelling (|Swinbank et all 
I2006D . 

There are two other observational links between these 
QSOs with M B « -25 (Mbh ~5x 1O 8 M ) and SMGs. 
Firstly, the QSOs have s pace densities at z ~ 2 some 5-10 
times lower than SMGs (|Bovle et al.l |2000| ; IChapman et al.l 
2003). This ratio of volume densities is roughly consistent 
with the ratio of the relative SMG and QSO lifetime, ~ 100— 
300 Myrs and ~ 20-40 Myrs, of ~ 5 times. Moreover, in P~3l 
we found a hint that this subset of our submm-detected 
QSOs lies on the present-day MBn/M sp h relation and simi- 
lar to SMGs, which is consistent with them being 'transition 
objects' between SMGs and QSOs in the proposed evolu- 
tionary sequence in which they would eventually evolve into 
passively evolving spheroids. 

Although this conclusion is based on a small sam- 
ple of QSOs, if the dynamical and gas properties of these 
three QSOs are representative of this subset of submm- 
detected QSOs, then our results are consistent with ~ L* , 
Mbh — 10 8 Mq submm-detected QSOs being 'transition 
objects' in the proposed evolutionary sequence, represent- 
ing a very brief prodigious star formation phase where the 
BH is undergoing rapid growth synchronously with the stel- 
lar mass. Submm-detected QSOs do not possess sufficiently 
large gas reservoirs (containing less gas than SMGs on av- 
erage) to sustain the implied SFRs for very long, providing 
a possible explanation for why some of our submm-detected 
QSOs are undetected in CO (a single weak CO detection 
and two non-detections of A/bh — 10 8 Mq submm-detected 
QSOs). 

Where do the optically luminous submm-detected 
QSOs fit into this evolutionary scheme? The average SMG 
would need ~ 1 Gyr (more than its lifetime) in order to 
grow their BHs by the required factor of £ 60 to host 
a Mbh — 4 x 10 9 M Q BH seen in the optically lumi- 
nous subset of submm-detected QSOs. In addition, these 
SMGs would need to begin with gas reservoirs containing 
Mgas Z 4 x 10 11 Mq in order to leave luminous QSOs with 
detectable reservoirs of residual gas. Alternatively, these rare 
optically-luminous QSOs could be periodically 'repeating' 
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SMGs, with SMG progenitors with M gas Z 10 11 M s and a 
surface density of ~ lOdeg -2 , which are consuming their 
gas and growing their BHs synchronously, but which accu- 
mulate additional BH mass and gas through merging with 
other S MGs in dense en vironments (i.e. damp 'dry merg- 
ing'; e.g. iBell et al"1l2004h . However, we currently have not 
observed any SMGs with such large gas masses, and given 
their rarity, we will require large scale submm blank-field 
surveys such as those planned with SCUBA-2 to find these 
potential progenitor SMGs. 

Hence our current results split the submm-detected 
QSOs in our survey into two categories (although a bigger 
sample would probably give a continuous trend): 

(1) Our observations of the optically luminous submm- 
detected QSOs indicate they have Mbh > 1O 9 M0. It is 
very difficult to grow the BHs in typical SMGs to such large 
masses in the typical lifetimes of SMGs or submm-detected 
QSOs. We therefore suggest that these QSOs cannot be re- 
lated to the evolutionary cycle of typical SMGs. Instead, 
these intrinsically rare QSOs most likely evolve from an 
equally rare subset of the SMG population - whose sur- 
face density is low enough that it has yet to be detected in 
current small-area submm surveys. Simply based upon the 
masses of their BHs, it is clear that these systems can only 
evolve into the most massive elliptical galaxies found in local 
rich clusters. 

(2) We have found a subset of our submm-detected 
QSOs with lower optical luminosities which have gas and 
BH masses and space densities which are consistent with 
them being potential 'transition objects' in the proposed 
evolutionary scenario where SMGs evolve into QSOs. This 
subset of submm-detected QSOs have Mbh — 10 8 Mq and 
gas consumption timescales consistent with a phase of prodi- 
gious star formation where the BH has grown substantially 
during the preceding ~ 200 Myr SMG phase. The fast build 
up of the stellar population is consistent with an evolu- 
tionary path ending in a population of present-day massive 
spher oids, where homogen eous old stellar populations are 
seen (|Swinbank et alj|2006h . 



6 CONCLUSIONS AND FUTURE WORK 

We have carried out a millimetre interferometry survey of 
nine submm-detected QSOs at z = 1.7-2.6 in order to test 
the link between SMGs and QSOs at the era where these two 
important populations were most numerous. We include in 
our analysis comparable observations of a similarly selected 
QSO from the literature to provide a final sample of ten 
submm-detected QSOs. To support this survey we obtained 
near-infrared spectroscopy of these QSOs to derive accurate 
systemic redshifts needed to tune the millimetre receivers 
to the correct frequencies. The near-infrared spectra also 
provide Ha fluxes and line widths needed to derive reliable 
BH mass estimates for the QSOs. Our main findings are: 

(1) We detect CO emission from six of the ten submm- 
detected QSOs in our sample, confirming that they con- 
tain a significant amount of molecular gas and that a large 
fraction of the mm emission is from starbursts. The me- 
dian gas mass of our sample (including non-detections) is 
(2.5 ± 0.7) x 10 10 M Q , similar to that found for z ~ 2 
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Figure 10. Af gas as a function of Mbh for our sample of CO- 
observed submm-detected QSOs. Recall that the quoted errors for 
M gas do not include the uncertainty on the CO-to-gas conversion 
factor. The arrows (offset slightly for clarity) indicate the move- 
ment of an average SMG in terms of its gas mass depletion and 
synchronous BH growth in arbitrary scalable timecales of 10 and 
20 Myr, assuming a SFR of 1000 Mq yr _1 and Eddington-limited 
BH growth. This demonstrates that an SMG could evolve into 
M BH 10 s M Q , M gas ~ 0.6 x 10 10 M Q submm-detected QSOs 
in a reasonable timescale, whereas an average SMG would need 
substantially (~ lOx) more gas and more time (Jt SMG lifetime) 
to evolve into the M B h - 4 X 10 9 Mq, M gas ~ 3.0 X 10 10 Mq 
submm-detected QSOs. This suggests that the submm-detected 
QSOs hosting BHs of M BH 10 8 M with M gas < 1 x 10 10 M Q 
comprise 'transition objects' that we can use to probe the inter- 
mediary evolutionary stage between the SMG and luminous QSO 
phases, while the rarer more luminous QSOs with Mbh ^ 10 9 Mq 
are not related to typical SMGs. 

SMGs and to z ^ 4 QSOs. The star formation efficien- 
cies of our QSOs are also comparable to those measured 
for SMGs, 250 ± 100 L (Kkms -1 pc 2 )" 1 , suggesting that 
the gross properties of the star formation in the QSOs are 
like those seen in SMGs. Adopting a 2kpc scale size for the 
gas distribution in the QSOs and a typical inclination of 
20° we derive a median dynamical mass of M(< 2kpc) ~ 
(2.1 ± 1.4) x 10 11 M , similar to SMGs (assuming an incli- 
nation angle appropriate for random inclinations). We find 
a lower incidence of double-peaked CO line profiles in the 
QSOs, compared to SMGs, which we believe results from a 
selection bias towards lower average inclination angles for 
the QSOs. 

(2) Our near-infrared spectroscopy indicates a median 
black hole mass in our QSO sample of (1.8 ± 1.3) x 10 9 M . 
Combined with our dynamical estimates of the spheroid 
mass, these yield M B H/M sph ~ 9 x 10~ 3 . This M B H/M ap h 
ratio for this sample of submm-detected QSOs at z = 2 is 
an order of magnitude larger than the local ratio, although 
M sp h suffers from large uncertainties due to the unknown 
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CO radii and inclination angles. This ratio is also signifi- 
cantly above that seen for SMGs at z ~ 2. However, this 
comparison masks a broad range in BH masses within our 
QSO sample and so we split the sample into two subsets 
based on their BH masses. 

(3) Looking at the optically luminous submm-detected 
QSOs in our sample we find that we detect CO emission in 
5/6 of these QSOs. However, the estimated BH masses for 
these QSOs, Mbh ^ 1O 9 -1O 1O M , are too large (and their 
number densities too small) for them to be related to typical 
SMGs in a simple evolutionary cycle. We propose that the 
progenitors of these most massive QSOs are a rare subset 
of SMGs with M gas > 4 x 10 11 M with a number density 
of ~ lOdeg -2 which will be possible to detect with future 
SCUBA-2 surveys. 

(4) For the optically less luminous (~ L*) submm- 
detected QSOs, we marginally detect one source in CO 
and obtain sensitive limits for three further QSOs. The BH 
masses for these systems are Mbh — 10 8 Mq , similar to the 
estimates for BHs in SMGs. These submm-detected QSOs 
are consistent with being 'transition' objects between SMGs 
and submm-undetected QSOs, as we show it is feasible to 
link their BH masses to those of SMGs by Eddington lim- 
ited growth for a period comparable to the gas depletion 
timescale of the QSOs, ~ 10 Myrs. The space density of these 
QSOs is also in rough agreement with that expected for the 
descendents of SMGs given current estimates of the relative 
lifetimes of QSOs and SMGs. We conclude that these ~ L* , 
Mbh > 10 s Mq submm-detected QSOs are consistent with 
being in a very brief prodigious star formation phase, and 
that they simply do not possess sufficiently large gas reser- 
voirs to sustain the SFR (which is why these might be less 
often detected in CO), although a larger sample of CO ob- 
servations of submm-detected QSOs with these BH masses 
is required for confirmation. 

To make further progress on understanding the evo- 
lutionary links between SMGs and QSOs requires a larger 
survey of the submm and CO emission from typical QSOs 
(Mb ~ —25). In addition, measurements of other CO transi- 
tions for the submm-detected QSOs (e.g. from IRAM 30-m, 
ALMA, EVLA, and SKA) are required to place better con- 
straints on the temperature and density of the molecular 
gas and thus provide a more accurate determination of the 
line luminosity ratios and hence total gas masses of these 
systems. Similarly, higher resolution CO observations are 
essential to put strong constraints on the reservoir sizes and 
inclination angles, and hence Md yn , needed to compare the 
two populations. Finally, better measurements of the far- 
infrared SEDs (with SABOCA, SCUBA-2 or Herschel) will 
yield more accurate measures of Lfir and Td us t for submm- 
detected QSOs to constrain the contribution from an AGN 
component. 
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